PROCEEDINGS 

°F SCIENCE 



o 
o 

(N 
> 

o 
oo 



> 

oo 
(N 
(N 



in 
o 



X 



Unifying the Fixed Order Evolution of Fragmentation 
Functions witli tlie l\/lodified Leading Logarithm 
Approximation 



Simon Albino"^ 

2nd Institute for Theoretical Physics, Hamburg University, Germany 
E-mail: 



simonSmail . desy . de 



Bernd Andreas Kniehl 

2nd Institute for Theoretical Physics, Hamburg University, Germany 

Gustav Kramer 

2nd Institute for Theoretical Physics, Hamburg University, Germany 

Wolfgang Ochs 

MPI fur Physik, Munich, Germany 
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the degrees of freedom given to the quark sector and the inclusion of the fixed order terms. We 
find that data from the largest x values to the peak region can be better fitted than with other 
approaches. 
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The current optimum description of single hadron inclusive production is provided by the QCD 
paiton model, which requires fragmentation functions (FFs) D^(;c, 2^) describing the probability 
for a parton a to emit a hadron h caiTying a fraction x of its momentum. Omitting the h and a labels, 
the evolution of the FFs in the factorization scale at large and intermediate x is well described 
[|] by the leading order (LO) DGLAP equation [g] 
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where P^^\z) are the LO splitting functions calculated from fixed order (FO) pQCD, and as{Q^) = 
oii(G^)/(2?r). As z — > 0, the LO splitting function asP^'-''\z) diverges due to terms of the form Us/z- 
These double logarithms (DLs) occur at all orders in the FO splitting function, being generally of 
the form {I / z){aslnz)^ {agln^ zY for r = — l,...,oo. As x decreases, Eq. (jl|) will therefore become 



a poor approximation once ln(l/x) = 0{a 
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A small X description is obtained by resumming 



these DLs using the Double Logarithmic Approximation (DLA) given by 
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for D = {PziDg), where Dz = -^.H'^/^iiDq + Dg) is the singlet FF. The evolution of the valence 
quark and non-singlet FFs vanishes in the DLA. 

We now construct an approach suitable for large and small x simply and consistently by using 
the DLA to resum the DLs in the DGLAP evolution, as described in more detail in Ref. [Q]. We 
will use Eq. ^ to modify a,P^'^^ in Eq. (0) to 

a,p(°)(z)^P°L(2,a,) + a,p(°)(2), (3) 

where P^^{z,as) contains the complete DL contribution, while asP^^\z) is obtained by subtracting 
the LO DLs, already accounted for in P^^, from asP^^\z) to prevent double counting. To obtain 
P^^, we work in Mellin space, /(ft)) = /(J dxx^f{x). Upon Mellin transformation, Eq. @ becomes 

(0 + 2^^]^^-2CAa,iQ^)A-lco + 2 
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a,(e2)p('')(«) 



D(ftj,e')=0, (4) 



where for completeness we have also accounted for a^P , which can be neglected in the following 
calculation of P^^. Making the replacement in Eq. (||) in Eq. (|l|), taking its Mellin transform. 



jj^D{(0, Q^) = P"^{(0,a,{Q^))D{(0, Q^), and then substituting this into Eq. @ gives 2{P"^y + 
^pDL _ iCaUsA = 0. We choose the solution P^^{co,a,) = j ( - (O + \/ (0^ + ISCAaA since its 



expansion in as yields at LO the result asP^^^'^^COjas) = 2CaA^, which agrees with the LO DLs 
from the literature [||]. The resummed result in x space is then P^^{z, as) = fl^^S^/^ (^4^CAasln jj , 
with Ji being the Bessel function of the first kind. 

If we approximate asP^^\(o) in Eq. by its SLs, defined at LO to be the coefficients of (O^, 

^Trhj, then if we apply the 
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approximate result that follows from the DLA at large Q, 
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the gluon component of Eq. becomes the MLLA differential equation [Q]. Therefore we con- 
clude that, since we do not use these two approximations, our approach is more complete and 
accurate than the MLLA. 

We now compare our approach to normalized differential cross section data for light charged 
hadron production from e^e^ {'Yi'Z') -^h + X. We fit the gluon Dg{x, QI) and the quark FFs 

Duc{x,Ql) = \ {Duix,Ql)+D,{x,Ql)) , Ddst{x,Ql) = ^ {Dd{x,Ql) + Ds{x,qI) + Di,{x,qI)) , 

(6) 

where Qo = 14 GeV. Since the hadron charge is summed over, we set Dq = Dq. For each of these 
three FFs, we choose the parameterization D{x,Q^) = A'^exp(— cln^;c);c"(l —x)^. To prevent too 
many free parameters, we use Eq. to fix c„f = Cdsb = Cg and a,,^ = OLdsb = otg- Since all scales 
are above the bottom quark mass, we set = 5. 

Performing a fit using the FO approach to LO, we obtain = 3 and the results in Fig. |l]. 
The fitted result of Aqcd = 388 MeV is quite consistent with that of other analyses, at least within 
the theoretical error of a factor of O ( 1 ) . It is cleai" that FO DGLAP evolution fails in the description 
of the peak region and shows a different trend outside the fit range. 




Figure 1: Fit to data in the FO approach to LO. Some of the data sets used for the fit are shown, together 
with their theoretical predictions from the results of the fit. Only data for which t, = ln(l /x) < In 
were used, indicated by the vertical lines. Each curve is shifted up by 0.8 for clarity. 

Now we perform the same fit again, but using our approach, i.e. Eq. (|l|) with the replacement 
in Eq. (^), for the evolution. The results are shown Fig. ^ We obtain Xbf = 2.1, a significant 
improvement to the fit above with FO DGLAP evolution. The data around the peak is now much 
better described. The energy dependence is well reproduced up to the largest ^/s value, = 202 
GeV. We obtain a rather large Aqcd = 801 MeV. We note that had we made the usual DLA (MLLA) 
choice Q = ^fsjl instead of our choice Q = ^/s as is done in analyses using the DGLAP equation, 
we would have obtained half this value for Aqcd- A treament to NLO is required to understand 
this problem, as well as treatment of hadron mass effects which are important at small x. 
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Figure 2: As in Fig. 1, but using our approach for the evolution. 



In conclusion, we have proposed a single unified scheme which can describe a larger range in 
X than either FO DGLAP evolution or the DLA. Our scheme allows a determination of quark and 
gluon FFs over a wider range of data than previously achieved, and should be incorporated into 
global fits of FFs such as that in Ref. [||] since the current range of 0. 1 < .x < 1 is very limited. 
Our approach should be expected to improve the description of other inclusive hadron production 
processes, e.g. those involving protons in the initial state. This work was supported in part by DFG 
through Grant No. KN 365/3-1 and by BMBF through Grant No. 05 HT4GUA/4. 
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